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Abstract. We present the new and updated BABAR measurements of CP-violation studies for 
many b — > s penguin decay modes. We report the first observation of mixing-induced CP-violation 

in B° — ^ rj K° with a significance (including systematic uncertainties) of 5.5<r. We also present 
the first observation of the decay B° — > p°K°. Using the time-dependent Dalitz plot analysis of 
B° — » K + K~ K° decay, the CP-parameters Acp and /3 e // are measured with 4.8<r significance, and 
we reject the solution near 7r/2 — /9 e // at 4.5a. We also present the updated measurements of CP- 
violating parameters for B° — > Kg-K , KgKgKg and -ir Tr Kg decays. An updated measurements 

of the CP- violating charge asymmetries for B — > r\ K ± , r]K ± , loK ± decays are also presented. 
The measurements are based on the data sample recorded at the T(4S) resonance with BABAR 
detector at the PEP-II B-meson Factory at SLAC. 
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1 Introduction 

Measurements of time-dependent CP-asymmetries in 
B° meson decays through Cabibbo-Kobayashi-Maskawa 
(CKM ) favoured b — ► ccs amplitudes [1 have provided 
crucial tests of the mechanism of CP- violation in the 
standard model (SM) [2]. A major goal of B Factory 
experiments is now to search for indirect evidence for 
New Physics (NP). One strategy is to compare the 
measured value of the CP-violation (CPV) parame- 
ters from b —>■ ccs to independent determinations of 
the same quantities, from processes that are sensitive 
to the contributions of NP effects through loop (pen- 
guin) diagrams. In the decay chain T(45) — > B°B , 
one of the two B mesons decays into a CP eigenstate 
fcp at time tcp, and other decays into a flavor spe- 
cific state /tag at time itag- The CPV parameters are 
measured using time-dependent decay rate 



x [1 + q{S f sm(Am d At) - C s cos(Am d At)}} (1) 

where At — tcp — itag, t is the B lifetime, Amd is 

the B°B° mixing frequency and q = +1(— 1) when 

f tag = B°(B ). The parameters Sf and Cf represent 
the mixing-induced and direct CPV respectively. For 
penguin-dominated processes b — > sqq (q — u,d,s), 
Sf and Cf are expected to be consistent with the 
values from b — > ccs decays ( namely, Cf ~ and 
Sf ~ sin 2(3). Additional CKM suppressed contribu- 
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tions to the amplitude can induce only small devia- 
tions from the expectation. On the other hand, addi- 
tional loop contributions from NP processes may pro- 
duce large observable deviations [3]. In this paper, we 
present new and updated measurements of b — > s pen- 
guin decays. 

2 Analysis Overview 

The updated and new measurements presented here 
are based on the data collected with the BABAR de- 
tector [J at the PEP-II asymmetric energy e + e~ col- 
lider, located at SLAC, USA. Unless explicitly stated, 
all the analyses are presented with a data sample cor- 
responds to 383±4 million BB pairs recorded at T(AS) 
resonance (Center-of-mass energy (CM) y/s = 10.58 
GeV) and 37 fb^ 1 data collected about 40 MeV below 
the T(45) resonance ( "off- resonance" ) . The BABAR 
detector provides the charged particle tracking through 
a combination of a 5-layer silicon vertex tracker (SVT) 
and a 40-layer drift chamber (DCH), both operating 
within a 1.5 T magnetic field generated by a supercon- 
ducting solenoidal magnet. Photons are identified in an 
electromagnetic calorimeter (EMC) surrounding a de- 
tector of internally reflected Cherenkov light (DIRC), 
which associates Cherenkov photons with tracks for 
particle identification (PID) . Muon candidates are iden- 
tified with the use of the instrumented flux return 
(IFR) of the solenoid. 

To discriminate between signal and background, 
two independent kinematic variables, the beam energy 
substituted mass toes = \/s/4 — (p*b) 2 an( i ener gy 
difference AE = E* B — y/s/2 are used. Here Eg and 
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p* B are the CM energy and 3-momentum of B candi- 
date, respectively. For signal decays, toes (AE) peak 
at nominal P-mass (at zero) with a resolution of a few 
(a few tens of) MeV. We use event topology to reject 
the continuum background (arises primarily from ran- 
dom combinations of particles in e + e~ — > qq events 
(q = u,d,s,c) events which have a "jet-like" struc- 
ture against the real B events which are more uniform. 
To improve the signal-background separation, some of 
these variables are combined in a multivariate algo- 
rithm (a Fisher discriminant or Neural Network etc.) 
which is trained on the signal monte carlo and qq back- 
ground and used in the fits. 

The physical quantities (e.g. signal yields or CP- 
parameters) are determined from extended and un- 
binncd maximum likelihood (ML) fits. The likelihood 
function include the probability density functions (PDF) 
whose shapes are based on MC or data (sidebands 
or off-resonance). Their most important parameters 
are left free. In addition to signal and continuum, B- 
background categories are also included in the fits. 
Then we calculate the systematic errors arising from 
several sources e.g. variation of signal PDF shape pa- 
rameters within their errors, P-background estima- 
tion, the effect of interference etc. In all the reported 
measurements of CP-violation below, the statistical 
errors dominates. 



iy + tt~K s contributions as a non-interfering background. 
The effects of interference with other resonances are es- 
timated and taken as a systematic uncertainties. The 
B° candidates are reconstructed from the combina- 
tions of p° — > ir + ir~ and Kg —> ir + ir~. An unbinned 
extended ML fit with the input of toes > AE, NN, cos 8 n + , 
TOtttt, At is used to extract the CP-asymmetry and 
branching fractions; here NN is the neural network to 
enhance discrimination between signal and continuum 
background, cos 9 n + is the angle between the Kg and 
the 7r + from the p° in the p° meson's CM frame and 
m n7r is the mass of pion pairs. The fit yields 111 ± 19 
signal events, we calculate the Branching Fractions 
(B° -> p°K°) = (4.9 ± 0.8 ± 0.9) x 10~ 6 , where the 
first error is statistical and the second systematic. The 
likelihood ratio between the fit result of 111 signal 
events and the null hypothesis of zero signal shows 
that this excluded at the 8. Oct level. When additive 
systematic effects are included, we exclude the null 
hypothesis at the 5.0a level. This is the first observa- 
tion of the P° p°K°. The fit for the CP-parameters 
gives S = 0.20±0.52±0.24 and C = 0.64±0.41±0.20. 
The detail of this analysis can be seen at Ref. [7]. 



3.3 B° -» K+K~K° 



3 Results 

3.1 B° -» rj'K 

Theoretically, it has the large Branching Fractions (65 x 
10~ 6 ) as compared to other penguin decay modes. 
Here we update our previous measurements. In addi- 
tion to the B° — > r\ Kg decays used previously, we now 
also include the decay B° — > r/ K\. The P-daughter 
candidates are reconstructed through their decays 7r° — > 
77, V -> 1l(Vu),V -> 7r + 7r~7r°(?737r), rf -> r/^ir+ir- , 
V ~ * ? y37r'7r + 7r ~, V — * p 7- Only 77 — > ?7 77 7r + 7r~ mode 
is used for the 77 K^. The detail of this analysis is 
described in Ref. [6]. We obtain the CPV parame- 
ters and signal yields for each channel from a ML 
fit with the input observables AE, toes, Fisher and 
At. Finally we obtain the CP-violations parameters, 
S = 0.58 ± 0.10 ± 0.03 and C = -0.16 ± 0.07 ± 0.03. 
The data used is 2.1 times as large as that of pre- 
vious measurements [5]. We observe mixing- induced 
CP- violation in B° — ► ij K° decay with a significance 
of 5.5 standard deviations (systematic uncertainties 
included). The results for direct CP- violation is 2.1 
standard deviations from zero. 

3.2 B° -> p°K° 

This analysis is reported with the 227 x 10 6 PP pairs. 
We take a quasi-two body approach, restricting our- 
selves to the region of the B° — > 7r + ir~K s Dalitz plot 
(DP) dominated by the p° and treating other P° — > 



Previous measurements of CP-asymmetries in B° — > 
K + K~ K° have been performed separately for the events 
K + K~ in the 6 mass region, and for events excluding 
the 4> region, neglecting interference effects among in- 
termediate states [8]. But now we have performed a 
time-dependent Dalitz plot analysis of B° — > K + K~K° 
decay from which we extract the values of CP- violation 
parameters by taking into account the complex ampli- 
tudes describing the entire B° and B Dalitz plots. We 
define the complex amplitude A for B° decay as a sum 
of isobar amplitudes [S] , 

A(m K + K ~, cos 8 h) = At 

r 

= J2 C r(l+brV { ^ +S r\f r ( mK+K -, C Os9 H ), (2) 

r 

where the parameters c r and <f> r are the magnitude 
and phase of the amplitude of the component r, and 
we allow for different isobar coefficients for B° de- 
cays through the asymmetry parameters b r and S r . 
Similarly, we can write the amplitude A for B° by 
replacing the plus signs with minus signs in above 
equation. Our isobar model includes resonant ampli- 
tudes </>, /o, Xco(lP), andX (1550) [10]; non-resonant 
terms; and incoherent terms of B° decay to D~K + 
and D~ K + . For each resonant term, the function f r = 
F r xT r x Z r describes the dynamical properties, where 
F r is the Blatt-Weisskkopf centrifugal barrier factor 
for the resonance decay vertex [11] , T r is the resonant 
mass line shape, and Z r describes the angular distribu- 
tion in the decay [12j . The CP- asymmetries parame- 
ters ^4cp,r(Rate Asymmetry) and P e ff,r (Phase asym- 
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Table 1. The fit fractions T T from the fit to the full B" 
K + K~K° Dalitz Plot. Errors are statistical only. 



Isobar Model 



tf>K° 
foK° 

X (1550)K° 

(K + K°K-) N r 

Xco(1P)K° 



12.5 ± 1.3 
40.2 ± 9.6 
4.1 ± 1.3 
112.0 ± 14.9 
3.0 ± 1.2 



Table 2. The CP-asymmetries for B° -> K+ K~ K° for 
the entire DP, in the high mass region, and for (f>K° and 
foK° in the low mass region. The first errors are statistical 
and second errors are systematics. 





A C p 


ft,/(rad.) 


Whole DP 


-0.015 ±0.077 ±0.053 


0.352 ± 0.076 ± 0.026 


High- mass 


-0.054 ±0.102 ±0.060 


0.436 ± 0.087+^^g 


4>K U 


-0.08 ±0.18 ±0.04 


0.11 ±0.14 ±0.06 


foK° 


0.41 ± 0.23 ± 0.07 


0.14 ±0.15 ±0.05 



metry) can be computed as 



.4 



CP,r 



\M 2 -\M 



-2b r 



\A r 



\A r \ 2 l + b^ 
Peff,r — + S r 



(3) 
(4) 



We first extract the values of the parameters of the am- 
plitude model, and measure the average CP-asymmetry 
in B° — > K + K~K° decay over the entire Dalitz plot. 
The fraction of events (T r ) obtained from the fit to 
whole Dalitz plot are shown in Table [TJ and the CP- 
asymmetries results are summarised in Table [2] The 
sum of the fractions is greater than one due to the 
interference. Using this model, we then measure the 
CP-asymmetries for 4>K° and foK° decay channels, 
from a "low mass" analysis of events with M K + K - < 
1.1 GcV/c 2 . Finally, we perform a "high mass" analy- 
sis to determine the average CP-asymmetry for events 
with M K + K - > 1.1 GeV/c 2 . The raw asymmetry be- 
tween B° - and B -tagged events in (a) the low-mass 
region and (b) high-mass region are shown in Fig- 
ure [l] and the obtained value of CP-parameters are 
given in Table[2] The values of (3 e ff and Acp obtained 
from the whole DP fit are consistent with SM expec- 
tations of j3 ~ 0.370 rad, Acp — 0. The significance of 
CP-violation is 4.8<j, and we reject the solution near 
7r/2 — P at 4.5cr. In the fit to the low-mass region of DP, 
we find f3 e f / lower than the SM expectation about by 
2a while in the high-mass region CP-parameters are 
compatible with SM expectations and we observe CP- 
violation at the level of 5.1a. The more detail of this 
analysis can be found in Ref. [13] . 




and -B°-tagged 



Fig. 1. The raw asymmetry between B° 
events in (a) low-mass region and (b) the high-mass region. 
The curves are projections of corresponding fit results. 



3.4 B° -» K s n° 

A precise measurement of direct CP-asymmetry and 
branching fraction for B° — > K°ir° is important as 
it is possible to combine direct CP-asymmetries and 
branching fraction measurements for all the four B — > 
Ktt modes to test precise sum rules. We have presented 
the updated measurement for this analysis. Since no 
charged particles are produced at the decay vertex of 
B so we determine the decay point by constraining 
the B production vertex to the interaction point (IP) 
in the plane orthogonal to the beam axis using only 
the Kg — > 7r + 7r~ trajectories. We account for possi- 
ble biases due to the vertexing technique by compar- 
ing fits to a large simulated sample of IP-constrained 
(neglecting the J /if) contribution to the vertex and us- 
ing Kg trajectory only) and nominal B° — ► J/ipKg 
events. The detailed analysis technique is discussed 
in Ref. [14]. We have measured the CPV parameters 
C K o„ = 0.24±0.15±0.03andS' K o 7r o = 0.40±0.23± 

0.03, and branching fraction B(B° -> K°ir°) = (10.34± 
0.66 ± 0.58) x 10 -6 . The first errors are statistical and 
the second are systematic. These values are consis- 
tent with SM predictions and the experimental value 
of sin 2/3; and the results supersede BABAR previous 
measurements. 



3.5 P° 



The decay B° — > KgKgKg is pure CP-even eigenstate 
and the theoretical uncertainty in the SM prediction 
of sin2/3 e // is less than 4%. The Ref. 112] gives the 
whole detail of this analysis. The whole data sample 
is divided into two subsamples, one where all three 

Kg mesons decays into the 7r + 7r~ channel (Pcp(h )) 

and another where one of the Kg mesons decays into 
the 7r°7r° channel, while the other two decay into the 
7r + 7r~ channel (Pcp(oo)) ■ We form ir° — * 77 candi- 
dates from pairs of photon candidates in the EMC. 

For each B C p(-\ ) candidate toes and AE are the 

two independent kinematical variable used while for 
B(jp(oo) reconstructed B° mass tub and the missing 
mass TO m i ss = y/ (q e + e - — qb) 2 is used. Here qh is 
the four-momentum of the Pcp(oo) candidate after a 
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mass constraint on the B meson has been applied. 
Since no charged particles are produced at the decay 
vertex of B, we have applied the same technique to 
find the vertex as used in B° — > Kgir (explained in 
Section 13. 4|) . The measured value of time-dependant 
CP-asymmetries, S = -0.71 ± 0.24 ± 0.04 and C = 
0.02 ± 0.21 ± 0.05, where the first errors are statistical 
and the second systematic. The statistical correlation 
between S and C is -14.1%. The statistical significance 
of CP-violation is 2.9a. These results agree well with 
the SM expectation and supersede our previous mea- 
surements. 

3.6 B° -> TrViTg 

The B° — * TT°TT°Kg final state is CP-even eigen state, 
regardless of any resonant substructure. Just like P° — ► 
K° s -k° and B° -> K° s KlK° s analyses, the same tech- 
nique is used to find the vertex of B candidates. We 
measure the CP- violating asymmetries in this decay 
from a sample of approximately 227 million BB pairs. 
From an unbinned extended ML fit we obtain S = 
0.72 ± 0.71 ± 0.08 and C = 0.23 ± 0.52 ± 0.13. When 
we fix the values of —S to the average sin 2/3 measured 
in b — > ccs modes, and C to zero, and re-fit the data 
sample the negative likelihood changes by 2.2ct. The 
detail of this analysis is described in Ref. [16] . 

3.7 B ± -> rj K^^K^^K* 

The decays B — > r\K are especially interesting since 
they are suppressed relative to the abundant B r\ K 
decays due to destructive interference between two 
penguin amplitudes. We present the updated measure- 
ment of CP- violating charge asymmetries (A c h) for 
B ± . For P ± — * decays, the direct CP- violation 
would correspond to a non-zero charge asymmetry de- 
fined as 

_ p(p- ^ ,n - p(p+ ^ /+) 

Two decay chains are reconstructed for the rj(— > 77 
and — ► 7T + 7r _ 7r°) and 77 (— > 7r + 7r~7r° andp°(—> 7r + 7r~)7) 
candidates. Then ML fit are done separately for each 
decay chain; the likelihood functions (including sys- 
tematic errors) are then combined to provide the charge 
asymmetries. We obtain the A c h(B + — > t]K + ) = — 0.22± 
0.11 ± 0.01, A ch (B+ -> 77'^+) = 0.010 ± 0.022 ± 0.006 
and Ach(B+ -> wK + ) = -0.01 ± 0.07 ± 0.01. We find 
no clear evidence for direct CP- violation charge asym- 
metries in these decays. The detail of this analysis can 
be found in Ref. [L7]. 

4 Conclusion 

We presented new and updated CP- violation stud- 
ies for many b — > s penguin decay modes. We re- 
ported the first observation of mixing-induced CP- 
violation in B° — > r\ K° with a significance (including 



systematic uncertainties) of 5.5cr. We also presented 
the first observation of the decay B° — ► p°K and a 
measurement of CP-violating parameters. Using time- 
dependent Dalitz plot analysis of B° — > K + K~K° de- 
cay, the Acp and /3 e // are measured. The significance 
of CP- violation is 4.8a, and we reject the solution near 
7r/2— /3 e f f at 4.5cr. We have also presented the updated 
measurements of CP-violating parameters for B° — » 
K%tt°, KlK° s K% and ir°ir°K% decays. All results su- 
persede previous BABAR results. Finally, an updated 
measurements of the CP- violating charge asymmetries 
for B^ — > 77 ,r\K^- ^wK^ decays are presented. The 
individual results of S ~ sin2/3 e y/ are consistent with 
world average of sin 2/3 from b — > ccs but naive aver- 
age of sin2/3 e // from all penguin modes is lower than 
sin 2/3 [18] which is opposite to theoretical prediction. 
Finally, we do not find any evidence of direct CP- 
violation in the decay modes reported. 
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